LUNFD6/(NFFL-7092) 1994 
hep-ph/9412323 



Rapidity-Dependent Low-]^^ Enhancement 



T. Csorgo 



1,2 



^MTA KFKIRMKI 
Q H-1 525 Budapest II4, P.O. Box 49. Hungary 

Q ' "^Department of Elementary Particle Physics, 



O 



X 



University of Lund, Sdlvegatan I4, S - 223 62 Lund, Sweden 



> : February 1, 2008 

m 
(N 
m 
(N 

0> : Abstract 

^ • The rapidity dependence of the \ow-pt enhancement is shown to be a sensitive 

. measure of the longitudinal source size for longitudinally expanding finite 



systems. 

Low-pi enhancement (LPTE) is conventionally understood as a rise in the ratio of the 
transverse momentum distribution of the particles produced in proton-nucleus {p + A) as 
well as heavy ion {A + B) reactions as compared to the transverse momentum distribution 
of the same particles in hadron-hadron reactions at the same center of mass energy per 
bombarding particle. An up-to-date summary of the status of the experimental data and 
theoretical understanding of LPTE can be found in ref. |1| which concluded: " In order 
to understand the results, theoretical models should ... explain the rapidity dependence 
observed in the data." 

In a recent paper the invariant momentum distribution (IMD) and the Bose-Einstein cor- 
relation function (BECF) has been calculated for longitudinally expanding finite systems . 



As an application of those results it is shown here that the longitudinal expansion together 
with the finite longitudinal size of the expanding tube naturally leads to a rapidity- dependent 
LPTE. 

Resonance decay effects and other, more exotic sources of LPTE, like the decay of 
droplets of Quark-Gluon Plasma or the effects of attractive potentials ^ leading to a 
change in the dispersion relation shall modify primarily the transverse momentum distribu- 
tion. However, these modifications are due to a mechanism which is independent (modulo 
isospin) of the longitudinal expansion thus the rapidity dependence of the LPTE does not 
naturally follows from these models. Partial thermal equilibrium due to a longitudinal ex- 
pansion was shown to lead to additional soft pions in ref. 0. These models explain neither 
the point why the enhancement is similar in p + A and A + A collisions nor the observed 
rapidity dependence of the LPTE data. 

Instead of studying the pt spectrum directly, let us concentrate on the coupling between 
the transverse momentum spectra and the longitudinal expansion. 

Longitudinally expanding finite systems. High energy heavy ion collisions may create 
one-dimensionally expanding systems, especially in the case of light projectiles. Heavier pro- 
jectiles may create three-dimensionally expanding systems. However, the three-dimensional 
expansion does not alter the coupling between transverse momentum spectra and rapid- 
ity distribution in the limit — m as long as the expansion remains in the same limit 
longitudinally boost-invariant. 

The four-momentum is given hj p = {E,p) = {E,px,Py,Pz)- The particle is on the mass 



shell = E"^ — p^, the transverse mass is denoted by rut = J E"^ — pI, the rapidity by 



y = 0.5 In 



^E+p,^ 
\E-PzJ 



the space-time rapidity by rj = 0.5 In 



\t-rzj 



and the transverse 



momentum is indicated by pt = ^p1 + Py. 

If a locally thermalized relativistic momentum distribution creates the correlations be- 
tween rapidity and space-time rapidity, the one-particle IMD can be calculated for arbitrary 
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space-time distribution functions as 

d Tl POO 

= gim) j_ dri G{ri) Jmt{v,y)- (1) 

dydrn^ °° 

In this equation the transverse and temporal components of the space-time distribution have 
already been integrated over p , and a finite size in the space-time rapidity is introduced by 
the distribution function G{r]). The various rapidity-independent modifications of the low-pt 
spectrum are taken into account through the factor g{mt) which shall be discussed later. 
Moreover, the local relativistic Bose-Einstein or Fermi-Dirac distributions for longitudinally 
expanding systems are given by 

/ da{r])-p 

JmAV,y)dv= • (2) 

(27r)^ exp {p-u{r])/T)±l 

Here / is the degeneracy factor, da{r]) ■ p is the inner-product of the volume-element of the 
freeze-out hypersurface and the four-momentum p = (mtCosh(?/), p^, Py, mjsinh(?/) ), while 
u{ri) stands for the four- velocity of the expanding matter on the freeze-out hypersurface and 
T is the freeze-out temperature. The + sign stands for fermions, the — sign for bosons. 
For the four-velocity of the ID expanding matter we assume the scaling Bjorken profile ^ 

as 

u{ri) = (cosh(?7), 0, 0, sinh(?7)). (3) 

Accordingly, Jrutiv^y) becomes a function of t] — y only. 

All the conventional and non-conventional modifications of the rrit distribution which 
are independent of the rapidity and aim at describing the low-pt behaviour at mid-rapidity, 
enter into the function g{mt). This function has infiuence on the rapidity dependence of the 
effective temperature. We assume that 

g{mt) = C exp{-mt/T^) (4) 

where T^, stands for a rapidity-independent eSective temperature in the low-pt limit. Such 
a structure is the simplest limiting distribution of many complex calculations attempting to 
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describe the LPTE. Note that here we are interested only in the \ow-pt behaviour and thus 
a deviation of the g{mt) distribution from the above simple exponential for medium or high 
values of rrit is not relevant for our considerations. Thus, the phenomenological parameter 
T^, stands for all the possible rapidity independent modifications of the rrit distribution at 
\ovf-pt, including decays of resonances, possible modifications of the dispersion relations and 
possible more exotic effects. 

In this model, the source at a given space-time - rapidity t] is emitting bosons with 
the mean rapidity (y) = rj. The width of the emission can be calculated by approximat- 
ing the Bose-Einstein and the Fermi-Dirac distributions with the relativistic Boltzmann- 
distribution, 

JmtiViV) ~ C'{mt) exp (-mt{y - r]f /T^ (5) 

where C'{mt) = C rrit ex.]i{—mt/T) is independent of rapidity and can thus be absorbed 
into the g{mt) factor. Only the leading order mtli] — y)^/T term has been kept in the Taylor 
expansion of the logarithm of J^f Terms of O {{rj — yY) are next to leading order since 
the Boltzmann approximation is valid for nit/T » 1 only. A numerical comparison of the 
approximated to its approximation indicates that the approximation is excellent for rrit > 3T 
and reasonable down to rrit > 1.2T. For certain classes of the models the freeze-out 

temperature T is smaller than the effective temperature at mid-rapidity, = T + Tq- 
Here the geometrical contribution to the effective temperature is given in the low-pt limit 
by Tg = mR^a? /tq where the finite transverse geometrical size is Rq-, the mean freeze-out 
time is denoted by tq and a stands for the flow gradinent in units of the mean freeze-out 



time. A preliminary analysis of the pion IMD [|rT|] indicates that vriT^jT ^ 1.5, for kaons 
mx/T ^ 6, thus the approximation (|^) is warranted. 

A Gaussian approximation is also applied for the distribution function of r/ 



1 

G{ri) = exp 



(6) 
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exp 



iv-y) 



2\ 



(7) 



where 



Ar]T (mt) = jT/mt, 



(8) 



the mid-rapidity is denoted by i/q and the width of the space-time rapidity distribution is 
given by Arj, which is a dimensionless measure of the longitudinal extension of the expanding 
system at freeze-out. By introducing this finite width Arj, we break the boost-invariance 
of our source in the longitudinal direction too. Thus we expect a non- stationary rapidity 
distribution. The flat rapidity distribution corresponds to the Arj oo limit. 
Gaussian results. The IMD can be calculated 



as 



d n 



dydm] {2-KAy^{mt)f'^ 



g{mt) exp( 



{y - y^f . 

2Ay\mt) 



(9) 



where 



Ay'^{mt) = At]^ + Ari^imt). 



(10) 



This relation introduces a coupling between the variables y and m^, i.e. their distributions 
do not factorize any more. At a fixed value of the rapidity, let us expand eq. (|^) around 
rrit = m as 



d n 



dy dm^ 



I 



(X exp 



nit 



V Uy) 



+ O {[nit - mf/rn^^ 



J 



which yields a Lorenzian effective temperature distribution in the low-pt region: 



Te{y) 



l + a{y- yof 



2m^ 



( 



Aif 



T 



m 



(11) 



(12) 



(13) 
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a dependence which is similar to the recent finding of the NA35 collaboration for the charged 
hadrons in S* + Pb collisions at CERN SPS [|T^. Such a change in the effective temperature 
with rapidity is a general property of heavy ion reactions performed both at CERN SPS 
and at the AGS, as measured by the NA35, E802/E859, E810, E814 collaborations valid 
pions, kaons protons and lambdas, see Fig. 3. in ref. |T3]. Thus the decrease of Te{y) in 
the target and projectile fragmentation region can be considered as a simple consequence 
of both the longitudinal expansion in the final stage and the rrit dependent width of the 
local thermal rapidity distribution. Simple kinematics results in a decrease of the average 
transverse momentum in the target and projectile fragmentation regions for p + p reactions 



too, see Fig. 9 of ref. [14 



With Te{y) given above the IMD can be integrated over rrit at a fixed y, resulting in 



dn Te{y) m + Te{y) 

— = Af exp 

dy T^, 171 + 



( 



{y - yof 



(14) 



\^ 2 [Ar]' + T/m) j 

where A/" is a normalization constant. This distribution describes for small values of a an 

approximately Gaussian rapidity distribution, with corrections which reduce the effective 

T 

width of the Gaussian to Ayl^^ < /S.rf H . 

m 

Eqs. ([T^JT^) are given in terms of three free parameters, such as the (dimensionless) 
longitudinal size Ai], the freeze-out temperature T and the effective temperature T* = Te{yo). 
Thus the longitudinal size At], which is hardly accessible to HBT measurements 0,0, can 
be determined from momentum space measurements alone! 

The experimental analysis of LPTE is based on a comparison of the IMD in p + A or 

A + B reactions to the IMD in p + p reaction at the same energy, 

-1 



R 



AB 
VP 



d nA+B 
dy dm\ 



d n, 



p+p 



\ 



^ dy dm\ j 



(15) 



This ratio can easily be evaluated from eq. (P) as 
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i?^/ = exp 



exp 



mt 



AB r-ppp 



X 



(y-yo)'(AV"-AV 



, 2AB\ 



\ 



(16) 



\^ (A^r^f^' + T/mt) {A^v^^ + T/rrit) J 
If T is independent of the reaction type, R^^ contains two factors. The first one is inde- 
pendent of the rapidity and may account for the various modifications of the pt spectrum at 
mid-rapidity due to rescattering, nuclear resonance production etc. This factor may yield 
the LPTE at mid-rapidity. The second factor determines the rapidity dependence of the 
effect. Generally, the degree of stopping in p reactions is smaller than in p -|- ^4 or ^4 -|- 5 
reactions at the same energy, due to the enhanced rescattering in the second and third case. 
This also implies Arj^ > Arj^"^ > Arj"^^, consequently indicate an increase in i?^/ within 
the target and projectile fragmentation region. Thus the LPTE in A + B or p + A reactions 
as compared to p -|- p reactions becomes a measure of the difference in the longitudinal sizes 
for the two reactions. 

Can the LPTE be utilized to measure the total longitudinal source sizes (not only differ- 
ences) ? 

Let us introduce as an auxiliary quantity the ratio of the IMD-s for rapidities yi and 1/2 

at the same nit the same reaction, either p + p 01 p + A or A + B 



R{yi,y2,mt) 



/ w2 \ 

a n 



^dy2dm^ J 



-1 



(17) 



dyi dm^. 

The rapidity-independent part of the momentum-distribution, g{mt) cancels from this ratio: 

driG{ri) Jmt{v,yi) 



R{yi,y2,mt) 



j dr]G{rj) J^, (77,^/2) 

Let us introduce the normalized ratio of the invariant momentum distributions. 



(18) 



LR(y,mt) 



(19) 



lim R{y,yo,mt) 



where yo stands, as before, for the mid-rapidity. This normahzed ratio shall be a measure 
of the LPTE as compared to the low-pt distribution at mid-rapidity. The normalized ratio 
Lpt has the property 

Jim^LR{y,mt) = LR{yo,mt) = 1. (20) 

Although the limit — > oo formally appears in the definition of L^, the limiting value 
of LR{y,mt) is 1, independently of the rapidity. This asymptotic value can be utilized for 
guiding the experimental normalization of this quantity. In the Gaussian approximation, 
Ln reads as 



+ 



LR{y,mt) = exp 



(21) 



The LPTE becomes by eqs. ([16| , |21|) exponentially enhanced in the target and projectile 



fragmentation region. This is due to the decrease ofTe{y) in the same region, cf. eq. (f2[). 
As a consequence, the low-p^ ratio Lfj{y,mt) becomes a sensitive measure of both the 
geometrical size of the emission region in the longitudinal direction Arj and the freeze-out 
temperature T. The LPTE has been calculated without explicit reference to the particle 
type, thus one might expect LPTE to exists for pions, kaons as well as for protons and any 
other heavier particles. At a given value of pt, heavy particles have a substancially reduced 
LPTE compared to pions according to eq. (|2TD . 

The expression ( ^I]) indicates that the rapidity-dependent LPTE is a consequence of the 
finite size of the longitudinally expanding system, because the enhancement vanishes in the 
limit the longitudinal size goes to infinity, 

lim Lii{y,mt) = l, (22) 

independently of y and nit. Thus it is more natural to compare the IMD, at a given rapidity, 
to the IMD at mid-rapidity in the same reaction, than to compare the IMD in p + A or A + B 
reactions to the IMD in p+p reaction, at the same energy. However, this latter possibility has 
been realized in most of the available data analysis on LPTE. As can be seen from eq. 
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the rapidity dependence of the ratio R^^ is also a consequence of the finite longitudinal 
sizes of the expanding systems since 



lim R^^ = exp 



vv 



( \ 

+ 



(23) 



which is rapidity independent. 

The results presented here are insensitive to resonance decay effects if these are not 
coupled to a specific rapidity region, since in such cases the rapidity independent effective 
temperature, may absorb their influence. A notable exception can be the delta reso- 
nance, which is copiously produced in the target and projectile fragmentation region. Since 
the effect of the longitudinal expansion is increasing with y — yo like an inverse Gaussian, 
eqs. ( |21|JT6| ), only a comparably strong coupling between the A production and the rapidity 



can influence it significantly. Such a strong coupling would be, however, rather unusual. 

The result for the LPTE is insensitive to the presence of the transverse flow, since similar 
results can be obtained for the IMD of ref. |T^, where the flow profile was approximated 



by a three-dimensional scaling fiow which describes a fully developed transverse fiow. The 
transverse fiow cancels since it infiuences the transverse momentum spectra only at high 
values of m^. 

It is straightforward to generalize the above results for arbitrary G{ri) functions, to be 
published elsewhere. 

Application. The E802 collaboration has found a characteristic bell-shaped curve for 
Te{y) in central S + Au reactions at 14.6 AGeV |T5|. It should be clear from the previous 



parts that the decrease of Te{y) with increasing \ y — yo \ results in a rapidity-dependent 
LPTE. Thus a fit to Te{y) is a kind of measure of the rapidity-dependece ot LPTE. The E802 
data are fitted with eqs. ([l2| , p^ in Fig. 1. The data are compatible with our description at a 



X^/NDF = 1.93. This result is somewhat surprising since we have utilized the scaling fiow 
profile to describe the longitudinal expansion. Note, however, that such a fiow profile may 
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develop not only in reactions where the incoming nuclei pass through each other, but also 
may develop in the final stage of the hydrodynamical evolution in the stopping region [jl6[ . 



The mid-rapidity at BNL AGS is at Umid ~ 1.7 for symmetric reactions, the fit yields 
somewhat smaller value of = 1-52 ± 0.25 for the asymmetric 5* + Au collision. The 
mid-rapidity temperature is found to be = 209 ±10 MeV, the freeze-out temperature 
cannot be determined precisely at the present level of the experimetal errors T = 60 ± 418 
MeV. Finally, we measure the length of the kaon-emitting volume at the freeze-out time, in 
dimensionless space-time rapidity units to be Ar^i^ = 0.39 ± 0.18. 

In summary, the model predicts a rapidity-dependent low-pj enhancement for p+p, p + A 
and A + B reactions, when compared to the transverse mass spectrum at mid-rapidity in 
the same reactions. The relative enhancement in p + A and A + B reactions compared to 
the transverse momentum spectrum in p + p reactions at the same energy turns out to be a 
consequence of the increased degree of stopping in p + A and A + B reactions as compared 
to the stopping in p + p collisions. 

The rapidity dependence of the low-pt enhancement can be used to measure the finite lon- 
gitudinal size for longitudinally expanding systems. This method is especially advantageous 
since the large longitudinal geometrical sizes of expanding systems appear in the radius pa- 
rameters of Bose-Einstein correlation measurements as correction terms only [@,^|lO|. Using 
a variation of this technique we have determined the length of the emitting volume at 
freeze-out in central S + Au reactions at AGS to be AriK+ = 0.39 ±0.18 space-time rapidity 
unit. 
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FIGURES 
S + Au at 14.6 AGeV, E802 data 
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FIG. 1. 



Fit to Te(y) of K'^ mesons in 14.6 AGeV 5 + Au reactions with eqs. (^JT^), solid line. Data 



points with error bars have been scanned from [15|. The fit parameters are Ar/ = 0.39 it 0.18, 
T, = 209 ± 10 MeV, T = 60 ± 418 MeV and yo = 1-52 ± 0.25 for the minimum of x^/NDF = 
5.8/(7 - 4) = 1.93 using tuk = 494 MeV. 
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